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Abstract. Let a be the scattering relation on a compact Riemannian manifold M with non- 
necessarily convex boundary, that maps initial points of geodesic rays on the boundary and initial 
directions to the outgoing point on the boundary and the outgoing direction. Let I be the length of 
that geodesic ray. We study the question of whether the metric g is uniquely determined, up to an 
isometry, by knowledge of a and I restricted on some subset T>. We allow possible conjugate points 
but we assume that the conormal bundle of the geodesies issued from T> covers T*M; and that 
those geodesies have no conjugate points. Under an additional topological assumption, we prove 
that a and £ restricted to T> uniquely recover an isometric copy of g locally near generic metrics, 
and in particular, near real analytic ones. 



1. Introduction and main results 

Let (M, g) be a compact Riemannian manifold with boundary. Let $* be the geodesic flow on 
TM, where for each (x,£) £ M, t i— > $ (x,£) is defined over its maximal interval containing t = 0, 
in particular this interval is allowed to be the zero point only. Let SM be the unit tangent bundle. 
Then dSM represents all elements in SM with a base point on dM. 

Denote 

(1) d+SM = {(x, € dSM; ±{v, £) < 0} , 

where v is the unit interior normal, (•, •) and stands for the inner product. The scattering relation 

(2) £ : dSM -» d+SM 

is defined by = (y, rj) = $ (x, £), where £ > is the first moment, at which the (unit speed) 

geodesic through (x, £) hits dM again. If such an C does not exist, we formally set £ = oo and we 
call the corresponding geodesic trapped. This defines also C(x, £) as a function C : d-SM — ► [0, oo]. 
Note that S and C are not necessarily continuous. 

It is convenient to think of E and C as defined on the whole dSM with S = Id and C = on 
d+SM. 

We parametrize the scattering relation in a way that makes it independent of pulling it back 
by diffeomorphisms fixing dM pointwise. Let k± : d±SM — > B(dM) be the orthogonal projection 
onto the (open) unit ball tangent bundle that extends continuously to the closure of d±SM. Then 
k± are homeomorphisms, and we set 

(3) a = k+ o S o kZ 1 : B(dM) — > B(dM), i = C o nZ 1 : B(dM) — ►[0,oo]. 

According to our convention, a = Id, I = on d{B{dM)) = S(dM). We equip B{dM) with the 
relative topology induced by T(dM), where neighborhoods of boundary points (those in S(dM)) 
are given by half-neighborhoods. 
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Let V be an open subset of B(dM). The lens rigidity question we study in this paper is the 
following: 

Given M and g\T(8M)> do a and I, restricted to T>, determine g uniquely, up to a pull back of a 
diffeomorphism that is identity on dM ? 

More generally, one can ask whether one can determine the topology of M as well. One motivation 
for the lens rigidity problem is the study of the inverse scattering problem for metric perturbations 
of the Laplacian. Suppose that we are in Euclidean space equipped with a Riemannian metric which 
is Euclidean outside a compact set. The inverse problem is to determine the Riemannian metric 
from the scattering operator, which is a Fourier integral operator, if the metric is non-trapping 
(see [Guj ). It was proven in |Guj that from the wave front set of the scattering operator, one can 
determine, under some conditions on the metric including non-trapping, the scattering relation on 
the boundary of a large ball. This uses high frequency information of the scattering operator. In the 
semiclassical setting, Alexandrova has shown that the scattering operator associated to potential 
and metric perturbations of the Euclidean Laplacian is a semiclassical Fourier integral operator 
that quantizes the scattering relation [Alj . [A2]. The scattering relation is also encoded in the 
hyperbolic Dirichlet to Neumann map on dM. Lens rigidity is also considered in [PoRj in the 
study of the AdS/CFT duality and holography, namely the idea that the "bulk" space-time can be 
captured by conformal field theory on a "holographic screen". The lens rigidity problem appears 
also naturally when considering rigidity questions in Riemannian geometry [CH IC2] . 

The lens rigidity problem is also closely related to the boundary rigidity problem. Denote by p g 
the distance function in the metric g. The boundary rigidity problem consists of whether p g (x,y), 
known for all x, y on dM, determines the metric uniquely. It is clear that any isometry which is the 
identity at the boundary will give rise to the same distance functions on the boundary. Therefore, 
the natural question is whether this is the only obstruction to unique identifiability of the metric. 
The boundary distance function only takes into account the shortest paths and it is easy to find 
counterexamples where p g does not carry any information about certain open subset of M, so one 
needs to pose some restrictions on the metric. One such condition is simplicity of the metric. 

Definition 1. We say that the Riemannian metric g is simple in M , if dM is strictly convex w.r.t. 
g, and for any x £ M , the exponential map exp^, : exp~ 1 (M) — > M is a diffeomorphism. 

The manifold (M, g) is called boundary rigid if one can determine the metric (and more generally, 
the topology) from the boundary distance function up to an isometry which is the identity at the 
boundary. It is a conjecture of Michel [Mi| that the simple manifolds are boundary rigid. This 
has been proved recently in two dimensions [PU] . for subdomains of Euclidean space [Grj or for 
metrics close to Euclidean [BIj . or symmetric spaces of negative curvature |BCG| . It was shown in 
[SU3] that metrics a priori close to a metric in a generic set, which includes real-analytic metrics, 
are boundary rigid. For other local results see [CDS] . [E] , [LSU] . [SU1| . The lens rigidity problem 
is equivalent to the boundary rigidity problem if the manifold is simple [Mi]. 

Of course there is more information in the lens rigidity problem if the manifold is not simple. 
Even so, the answer to the lens rigidity problem, even when T> = B(dM), is negative, as shown 
by the examples in [CK]. Note that in these examples the manifold is trapping, that is, there are 
geodesies of infinite length. The natural conjecture is that the scattering relation for non-trapping 
manifolds determines the metric uniquely up to the natural obstruction ([U]). There are very few 
results about this problem when the manifold is not simple. Croke has shown that if a manifold is 
lens rigid, a finite quotient of it is also lens rigid |C2| . 

This kind of data is overdetermined. At least when (M, g) is simple, knowledge of x, y, I (the 
graph of the boundary distance function) determines £, rj uniquely [Mi]. This can be extended, 
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under some non-degenerate assumptions, to the general case, see section[3j Instead of the boundary 
distance function, we need to know the multiple travel times between boundary points. Also, for 
simple manifolds, a known on the whole d_SM determines i uniquely. 

We redefine a in a way that removes the need to know g on T(dM). Denote by T°(dM) the 
tangent bundle T{dM) considered as a conic set, i.e., vectors with the same direction in T(dM) are 
identified. For any metric g\T{8M)i T°(dM) \0 is isomorphic to the unit tangent bundle S(dM) (in 
the metric g) but has the advantage to be independent of the choice of g. Given ^ f 6 B x (dM), 
we set 



i.e., A and 9 are polar coordinates of If £' = 0, then 9 is undefined. If £' = re±(£), knowing 
A, 9 is equivalent to knowing the angle that £ makes with the boundary, and the direction of the 
tangential projection of the same vector. Given two metrics g and g on M, and (a?, £') £ B(dM), 

0,O € B(dM), where B(dM) is related to g, we say that £' = |' iff = and £' = s|' for 
some s > 0. In other words, we require that £' and £' have the same polar coordinates (HJ). Note 

that this induces a homeomorphism B{dM) \— > B(dM) given by i— > \£'\g£' /\£'\g if £' 7^ 0, i— > 0. 
With that identification of B{dM) for different metrics, it makes sense to study a restricted to the 
same set T> for a family of metrics, and in particular, an a priori knowledge of g on T(dM) is not 
needed to define V and a on it. If a(x, £') = (y, r/), we just think of £' and 7/ as expressed in the 
polar coordinates (fjl|). Also, the notion of T> being open is independent of g. 

A linearization of the boundary rigidity problem and the lens rigidity problem, see section 14.31 
is the following integral geometry problem. Given a family of geodesies V with endpoints on dM, 
we define the ray transform 



of symmetric 2-tensor fields / (playing the role of the variation of the metric g), where (/, 9 2 ) is the 
action of / on the vector 9. Locally, (f,9 2 ) = fij9 l 9^. We will omit the subscript T to denote an 
integral over a chosen geodesic, or over all geodesies. Any such / can be decomposed orthogonally 
into a potential part dv and a solenoidal one f s (see section U]), and / vanishes on potential tensors. 
The linearized boundary rigidity or lens rigidity problem then is the following: can we recover 
uniquely the solenoidal projection f s of / from its ray transform? If so, we call ip s-injective. 

S-injectivity of I was proved in |PS] for metrics with negative curvature, in [Shl[ ISh2l |D| IPe] 
for metrics with small curvature and in [ShU] for Riemannian surfaces with no focal points. A 
conditional and non-sharp stability estimate for metrics with small curvature is also established in 
[Shi]. In [SU2] . we proved stability estimates for s-injective simple metrics, see ([26]) ; and sharp 
estimates about the recovery of a 1-form / = fjdx 3 and a function / from the associated If. These 
stability estimates were used in [SU2j to prove local uniqueness for the boundary rigidity problem 
near any simple metric g with s-injective /. In |SU3| . we showed that the simple metrics g for which 
/ is s-injective is generic, and applied this to the boundary rigidity problem. We note that in all 
the above mentioned results the metric has no conjugate points. On the other hand, in [SU4J we 
proved generic s-injectivity for a class of non-simple manifolds described below. Within that class, 
the boundary is not necessarily strictly convex, we might have conjugate points on the metric, the 
manifold might be trapping, and we have partial or incomplete information, i.e., we do not know 
the scattering relation or the ray transform for all sets of geodesies. 



(4) 



A=ie'i 9 e[o,i], 9 = e/\e\ a cT°(dM), 



(5) 
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Given (x, £) S V, let 7„-i/ x ^ denote the geodesic issued from «_ (x, £) with endpoint 7r(o"(x, £)), 
where 7r is the natural projection onto the base point. With some abuse of notation, we define 



Definition 2. We say that T> is complete for the metric g, if for any (z, £) E T*M there exists a 
maximal in M, finite length unit speed geodesic 7 : [0, 1] — > M through z, normal to such that 



We call the C k metric g regular, if a complete set T> exists, i.e., if B(dM) is complete. 

If z € dM and £ is conormal to dM, then 7 may reduce to one point. Since (JBJ) includes 
points where 7 is tangent to dM, and a = Id, £ = there, knowing a and I on them provides 
no information about the metric g. On the other hand, we require below that T> is open, so the 
purpose of © is to make sure that we know a near such tangent points. 

Definition 3. We say that (M,g) satisfies the Topological Condition (T) if any path in M 
connecting two boundary points is homotopic to a polygon c\ U 71 U ci U 72 U • • • U jk U ct+i with the 
properties that for any j, 

(i) Cj is a path on dM; 

(ii) 7j : [0,lj] — > M is a geodesic lying in M mt with the exception of its endpoints and is 
transversal to dM at both ends; moreover, k_ (t?(0), 77 (0)) € T>; 

Notice that (T) is an open condition w.r.t. g, i.e., it is preserved under small C 2 perturbations 
of g. 

We showed in [SU4j that if T> is complete, then Ijyf recovers the singularities of f s . Next, see 
also Theorem [2] below, under the same conditions, and assuming (T) as well, Id is s-injective for 
real-analytic metrics, and if k S> 2, also for generic metrics. 

To define the C k (M) norm in a unique way, and to make sense of real analytic g } s, we choose 
and fix a finite real analytic atlas on M. 

Theorem[T]below says, loosely speaking, that for the classes of manifolds and metrics we study, the 
uniqueness question for the non-linear lens rigidity problem can be answered locally by linearization. 
This is a non-trivial implicit function type of theorem however because our success heavily depends 
on the a priori stability estimate that the s-injectivity of Id implies, and the latter is based on 
certain hypoelliptic properties of Id, as shown in |SU4j . see (|26p . We work with two metrics g 
and g; and will denote objects related to g by a, £, etc. Note that (T) is not assumed in the first 
theorem. 

Theorem 1. Let go € C k {M) be a regular Riemannian metric on M with k ^> 2 depending on 
dim(M) only. Let T> be open and complete for go, and assume that there exists D' <s= T> so that 
Ig ,V * s s-injective. Then there exists e > 0, such that for any two metrics g, g satisfying 

( 8 ) \\g - 9o\\c k {M) + \\g - 9o\\c k (M) < e i 



Id(x,0 = I(% 



(6) 
(7) 





£ 



M fixing the boundary such that 



on V 



9 = ^*g- 
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Next theorem is a version of |SU4[ Theorem 3]. It states that the requirement that I g0t x>' is 
s-injective is a generic one for go. 

Theorem 2. Let Q C C k (M), with k S> 2 depending on dim(M) only, be an open set of regular 
Riemannian metrics on M such that (T) is satisfied for each one of them. Let the set V C B{dM) 
be open and complete for each g G Q . Then there exists an open and dense subset Q s of Q such that 
Ig,W * s s-injective for any g G Q s . 

Theorems Q] and [2] combined imply that there is local uniqueness, up to isometry, near a generic 
set of regular metrics. 

Corollary 1. Let T>' , Q, Q s be as in Theorem^ and let T> 3) T>' be open and complete for any 
g € Q. Then the conclusion of Theorem^ holds for any go G Q s . 

In section[3]we also prove that one can recover the jet of g on dM in boundary normal coordinates 
from a, I under a non-conjugacy assumption. There is no generic assumption, and the recovery is 
actually explicit. 

Remark 1 . Condition (T) in Theorem [21 and Corollary [T] in some cases can be replaced by the 
assumption that (M,g) can be extended to (M,g) that satisfies (T). One such case is if (M,g) is a 
simple manifold, and we study a, t on its maximal domain, i.e., T> = B{dM). In particular, we get 
local generic lens rigidity for subdomains of simple manifolds when T> is maximal. See section [5] for 
more details. 

2. Preliminaries 

We allow the geodesies to have segments on dM, then they are called geodesies if they satisfy 
the geodesic equation in a half-neighborhood of each boundary point. Such segments are included 
in determining the maximal interval, where $*(;c,£) is defined. If (M,g) is any extension of (M,g) 
(a Riemannian manifold of the same dimension of which M is a submanifold) , then any geodesic 
in M restricted to M is a geodesic in M. For the maximal geodesies in M we have the following 
property that indicates that the property of 7 to be maximal in M does not change under such 
extensions. 

Lemma 1. Let j(t), 0<t<l,0<l< 00, be a maximal geodesic in M. Let (M,g) be any C 1 ' 
extension of (M,g). Then there is no interval I D [0,1] strictly larger than [0,1], such that 7 can 
be extended as a geodesic in M for t G I, and {^(t); t G 1} C M. 

Proof. Suppose that there is such /. Without loss of generality, we may assume that J D [0, 1 + S], 
5 > 0. Then 7 solves the geodesic equation for t G [0,/ + S], with the Christoffel symbols T^- 
depending on g and its first derivatives restricted to M. Since g is continuous and has continuous 
first derivatives across dM, their restriction to M depends on g only. Therefore, 7 is a geodesic in 
M for t G I, and this is a contradiction. □ 

In [SU4] we studied geodesies originating from points outside M given some extension (M,g). 
We connect the notion of T> being open with the analysis in [SU4] by the following. 

Lemma 2. Let M be an extension of M, and let g be a C 1 ' extension of g on M. Let 70 : [0, 1] 1— > 
M be a unit speed geodesic with endpoints in M mt \ M such that 

(9) {( 7o (i),7o(t)); te[o,i]}ndISM CkZ^V). 

Then there exists a neighborhood W of (xo,£o) = (70(0), 70 (0)) such that any geodesic 7 with initial 
conditions in W and the same interval of definition t G [0,/] still satisfies Q. 
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Moreover, if g belongs to a class of extensions satisfying \\g\\(ji,i < A with some A > 0, then W 
can be chosen independently of g. 

Proof. Let Eq be the l.h.s. of ([9]). Then Eq is compact in SM. Given any neighborhood U of Eq (in 
SM), the set dSM \ U is compact, therefore there exists a neighborhood W of (xq,^q) such that 
the geodesic flow through W for < t < I will miss that set, therefore, its only common points with 
d-SM must be in U. To construct U, we first choose a neighborhood U X £ in SM of each point 
(x,0 G kZ X (V) so that U X £ ndlSM C k! 1 ^). This is easy to do in local coordinates because T> 
is open. Then we set U = U, £) gK -i(£>)^x,£- 

To prove the second part, we use the theorem of continuous dependence of solutions of an ODE, 
over a fixed interval, on the initial conditions and on the coefficients of the ODE. As long the 
Lipschitz constant related to the generator of the geodesic flow is uniformly bounded, we can 
choose W uniformly w.r.t. g. □ 

Lemma 3. Let (M,g), (M,g), 70 be as in Lemma\^ Let H be a hypersurface through xq = 70(0) 
transversal to 70, and set 7i = SM n 7r _1 (//). Then there exists a small enough neighborhood U 
of (xoj^o) i n *H> so that the geodesies issued from U with interval of definition t G [0,/] satisfy (0|) 
and are transversal to dM at all common points with it except for a closed set of initial conditions 
in U of measure zero. 

Proof. We only need to prove the statement about the measure zero set. We study the geodesic 
flow on SM issued from U. The corresponding geodesic is tangent to dM at some point x if the 
corresponding integral curve in the phase space is tangent to dSM at tt~ 1 (x). By Sard's theorem, 
this happens only on a closed set of measure zero. □ 

Those lemmas will be used to reformulate the results in |SU4j in the situation in the paper. 
In [SU4j . we extended the geodesies slightly outside M and parametrized them by initial points 
and directions on surfaces H m transversal to them, see section 14. 2\ instead by points and direc- 
tions on dM. This can be done, because when studying ir> the metric g is known and can be 
extended outside M in a known way. The reason for doing this was to prevent working with a 
parametrization, where the geodesies can be tangent to the surface, which is the case with dM. We 
can still do this even for the lens rigidity problem, using the boundary recovery result in section [3] 
below, see Proposition [2l We prefer however to parametrize the scattering relation by points on 
dM and corresponding directions. This does not preserve the smooth structure of the previous 
parametrization but the lemmas above show that it preserves the topology, at least. 

3. Recovery of the jet of g on dM. 

If dM is convex in a neighborhood of some point xq, then it is known that the jet of g near xq in 
boundary normal coordinates is uniquely determined by studying the "short" geodesies connecting 
xo and y, and then letting y — > xo, see [LSU| . This argument does not provide explicit recovery 
however. If dM is strictly convex near xq, then there is also a conditional Lipschitz stability 
estimate, see [SU3| . Below we show that we can do this also without the convexity assumption 
by studying "long geodesies" that do not converge to a point. A non-conjugacy condition will 
be imposed. Note that there is no generic assumption in the theorem below, and the proof gives 
explicit recovery. 

Let (M,g) be a smooth extension of (M,g). Let xo £ M, yo G M be endpoints of a geodesic 
70 : [0, 1] — >G M, and assume that xo, yo, are not conjugate points on 70. Then there exist 
neighborhoods U 3 xo, V 3 yo in M, such that for x G U, the exponential map neigh(7o(0)) 9 £ 1— > 
y = exp x £ G V is a local diffeomorphism, therefore it has a smooth inverse V 3 y 1— ► £ := exp" 1 y G 
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T X M such that exp" 1 yo = 70 (0). The geodesic 7^ : [0,1] — > M then connects x and y and is 
unique among the geodesies issued from x in directions close enough to £. One can also define a 
smooth travel time function t(x, y) on U x V by r(x,y) = [exp" 1 ?/^ with r(a;o,yo) = j'V(O) | ff . If 
there are no conjugate points on 70, then t locally minimizes the distance but in general, this is 
no longer true. On the other hand, r is always a critical value of the length functional and of the 
energy functional. 

In the situation described in the paragraph above, we will call y £ V visible from x 6 U, if 
7 a!) £ C M. In the next theorem, given G TM, we call y E M reachable from if there 

exists s > 0, such that 7a;,f (s) € M for f G [0, s], and 7a;^(s) = y. 

The so defined function r solves the eikonal equation (grad^, = 1 in U despite the possible 
existence of pairs of conjugate points not in U x V. Indeed, clearly, dr('y x ^(t),y)/dt\t=o = — 
Next, by the Gauss lemma, the x-derivatives of r(x,y) in directions perpendicular to £ vanish. 
Therefore, (£, grader) = — \£\ g , and the two vectors are parallel, therefore, grader = — and 
Igrad^ r\g = 1. Also, if 77 = grader, then for some I > 0, <l? z (x,£) = (y, r/). In particular, 7/ = grader, 
where the prime stands for a tangential projection. If <r is locally known, then integrating the known 
rj' = rf{y) along a curve on dM connecting y and yo recovers uniquely r(xo,y) up to a constant. 

Assume further that 70 is transversal to dM at both ends and does not touch dM elsewhere. 
Then 

(10) a(x,-gr&d' x T(x,y)) = (y,gr&d! y (x,y)), £(x, -grad^.r(x, y)) =r(x,y), 

Therefore, at least in this non-degenerate situation, knowledge of r(x,y) recovers uniquely a, i 
locally. 

The travel (arrival) times are widely used in the applied literature. Assume that g on T{dM) 
is known and fixed. Fix x, y on dM. We call the number r > a travel time between x and y, if 
there exists a geodesic of length r connecting x and y. Then we have a map that associates to any 
(x, y) € dM x dM a subset of (0, 00]. In the situation above, r(x,y) is one of the possible travel 
times between x and y. 

Theorem 3. Let (M,g) be a compact Riemannian manifold with boundary. Let (xo,£o) ^ S{dM) 
be such that the maximal geodesic 70 through it is of finite length, and assume that xq is not 
conjugate to any point in 70 Pi dM. If a and £ are known on some neighborhood of (xo,£o)> then 
the jet of g at xq in boundary normal coordinates is determined uniquely. 

Proof. To make the arguments below more transparent, assume that the geodesic 70 issued from 
0eO)£o) hits dM for the first time transversally at 70(^0) = yo> ^0 > 0. If 1$ = 0, the results is known 
[LSU], as we pointed out above (but this proof applies to this case as well). Then yo is the only point 
on dM reachable from (xo>£o)> and xq, yo are not conjugate points on 70 by assumption. Assume 
also that dM is strictly convex at xq, when viewed from the interior, i.e., the second fundamental 
form at xq is strictly negative. Then there is a half neighborhood V of xq on dM visible from yo 
[Sh2j. The latter is not always true for example when 70 is tangent to dM of infinite order at xq. 

Choose local boundary normal coordinates near xo and yo, and let go be the Euclidean metric in 
each of them w.r.t. to the so chosen coordinates. In boundary normal coordinates, o~, £ determine 
uniquely S, C. We can then consider a representation of E, denoted by E" below, defined locally 
on R n_1 x S n ~ l , with values on another copy of the same space. If (x,9) € R n_1 x S" 1 ^ 1 , then the 
associated vector at x G dM is £ = 0/\9\ g ; and T,*(x,9) = E(x,£). The same applies to the second 
component of E"(x, 9). Namely, if (y, if) = E(x, £), then we set u> = T)/\rf\ go , then E^ : (x, 9) \— * (y, u). 
Similarly, we set &{x, 9) = C(x, £). Let also #0 and lvq correspond to £0 and 770, respectively, where 
E(x ,£o) = (yo,Vo)- 
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Set t(x) := t(x,i/q), where r is the travel time function described above defined in a small 
neighborhood U of Xo in M. In the normal boundary coordinates x = (x', x n ) near xo, gin = $in, Vi. 
Since xo and yo are not conjugate, for r\ G S yo M close enough to rjo, the map r] i— > x € <9M is a local 
diffeomorphism as long as the geodesic connecting x and yo is n °t tangent to dM at x. Moreover, 
that map is known because it is determined by the inverse of £ near (xo,£o)- Similarly, the map 
S"™ -1 B uj i— > x is a local diffeomorphism and is also known. Then we know (x, —6) = Y,^(yo, —uj), 
and we know £^(yo, —u) = C^(x,9) = r(x). Then we can recover grad'r = —0'/\9\ g . Taking the 
limit uj — > u!q, we recover \9q\^ = g a p0Q6Q. We use now the fact that a symmetric n x n tensor fij 
can be recovered by knowledge of fijvtvl for N = n(n + l)/2 "generic" vectors k = 1, ... , iV; 
and such iV vectors exist in any open set on S"™ -1 , see e.g. |SU4| . Thus choosing appropriate 
n(n — l)/2 perturbations of 9qs, we recover g(xo). Thus, we recover g in a neighborhood of xq as 
well; we can assume that V is covered by that neighborhood. 

Note that we know now all tangential derivatives of g in V B xq. We showed above that r solves 
the eikonal equation 

(11) g al3 T xa T x e + Tin = 1. 

Next, in V, we know T x a, a < n — 1, we know g, therefore by we g is easy to 

see that r^n < on the visible part V, so we recover r^n there. We therefore know the tangential 
derivatives of r x n on dM near xo- 

Differentiate (fTTI) w.r.t. x n at x = xq to get 



(12) 



dg a P 



= 0. 

=.'■() 



Since 70 is tangent to dM at xq, we have r x n(xo) = by (fTT|) . The third term in the l.h.s. of (fT2j) 
therefore vanishes. Therefore the only unknown term in (|12p is G a @ := dg a " /dx n at x = xo- Since 
^(^o) = — £o> using the fact that gradr(xo) = — Co again, we get that we have to determine G al3 
from G a /3^Q^Q. This is possible if we repeat the construction and replace £0 by a finite number of 
vectors near £q, as above. So we get an explicit recovery of dg/dx n \oM in fact. 

Next, for x E V but not on dV, we can recover T x n x n(x) by (fT2|) because tv^x) < 0. By 
continuity, we recover T x n x n(xo), therefore we know r x n x n near xq, and all tangential derivatives of 
the latter. 

We differentiate (1121) w.r.t. x n again, and as above, recover d 2 g/d(x n ) 2 \dM near xq. Then we 
recover 9 3 r/0(x n ) 3 , etc. 

Let us return to the general case. We will show first that the assumptions about (xq,£o) are 
preserved under a small perturbation of that point. Let U be a neighborhood of (xo,£o) i n TM. 
Since 70 D dM consists of points that are not conjugate to xo, there is an open set W D 70 n dM 
in M that stays away from the points conjugate to x along the geodesies t 1— ► exp^ t£, < t < Iq if 
(x,£) € U and if U is small enough. Here Iq is the length of 70. On the other hand, the possible 
common points of those geodesies must be in W, if U is small enough, as in the proof of Lemma [2j 
By Lemma [TJ there exists s > 1 close enough to 1 so that (xo,s^o) £ U, and exp Xo (s^o) ^ M. 
The geodesies issued from (x, £) € ?7 close enough to (xo,s^o), and the same time interval [0, Iq] of 
definition, still have endpoints outside M, therefore they are longer than the maximal segment in 
M. On the other hand, we showed that they meet dM at points that are not conjugate to xo- We 
can now replace s£o by £0 and rescale the corresponding geodesic. This shows that the assumption 
in the theorem is preserved under a small enough perturbation of (xo,£o)- 

Let £ E = £0 + £V i < £ *S 1, where v is the interior unit normal at xq. In the coordinate system 
above, v = (0, . . . ,0, 1). Let j e be the geodesic issued from (xo,£e) until it hits dM for the first 
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time. By a compactness argument, we can choose a sequence Sj \ such that y £ . — » j/q that can 
be different from y but we still have that y$ € dM and ?/q is reachable from (xo, £o)- Then Z £j . — > 
with some Zq > 0. By assumption, xq and t/q are not conjugate points on 70, and the function 
t(x,u) is then well-defined near (xq,i/q) satisfying the eikonal equation. Next, r y £j connects xq 
and y £j , and only the endpoints are not in M mt . The advantage now is that ■j £j hits dM at xq 
transversely but it might be tangent to dM at y e . . 

For a fixed j, by [Sh2j (see the proof of Lemma 2.3 there) if U is a small enough neighborhood of 
xq on dM, then U can be expressed as the disjoint union U + U H U where -ff is a hypersurface 
on dM through xo, and at least one of the half-neighborhoods C/ ± is visible from y £j , let us say 
that this is U + . If 7 £j is transversal to dM at y ej , then even better, the whole U is visible from y £j 
if U is small enough. If n = 2, then if reduces to the point xq and the modifications are obvious. 
We set t(x) = r(x,y £j ). By the arguments in the special case above, we know t(x) on dM. 

We first recover g on dM in boundary normal coordinates. Let 9j = £ £j /|£ £j | So , and #0 = 
£o/|£o|g - As m t ne beginning of the proof, we take a one-sided derivative at xo, from U + , to 
recover —9 ej /\9 ej \g = — 9'o/\9o\ g + 0(sj). Take the limit £j — > to recover |#o|g- Thus we get g at 
xo and therefore, near xq. 

As above, from the eikonal equation, satisfied in U + , where r is defined and known up to a 
constant, if U + is small enough, we recover r^n on U + , and therefore rfn(xo) by continuity. Since 
T~x n (%o) = ~ £ j/\Cej \g, and therefore, r x « < near xo, we recover r x n(xo), and therefore, the normal 
derivative of r on dM on C/ + U -ff, and the tangential derivatives of the latter. 

By m, 

G a p(x )& $ /\£ £j \ 2 = (-2g af3 T xax nT x p + 2T x n. x n£j/\£ £ \g) 

\ / X=XQ 

where G aP = d a/3 /dx n . Since r is smooth near (xo,yg)' then T x n x n(xo) = T x ^ x n(xo,y £j ) remains 
bounded and even has a limit, as j — > 00. Similarly the other terms above have a limit. Therefore, 
the second term on the r.h.s. above tends to zero, as j — > 00, and we recover G a p(xo)^Q £q . 

To recover dg/dx n at xo, we need to perturb £q- As we showed above, (xq,C) satisfies the 
assumptions of the theorem, if £ is close enough to £o- We can therefore recover G a ^(xo)^ a ^ 13 for 
£ in a neighborhood of £oj which recovers G Q/ g(xo), and therefore dg/dx n at x = xo- 

To recover the higher order derivatives, we proceed in the same way. To recover d k g/d(x n ) k , 
we differentiate (jlip k times, and solve for d k g a @ / d(x n ) k at x = xq- The only unknown term in 
the r.h.s. will be d k+1 T / d(x n ) k+l but it will be multiplied by r x n that equals Sj/\^ £j \ g at xo- Then 
taking the limit j — > 00 will recover d k g/d(x n ) k as above. Then we recover d k+1 r / d(x n ) k+l needed 
for the next step, etc. □ 

Remark 2. If g has a finite smoothness g G C k {M), then the proof above implies that we can 
recover (explicitly) d a g\oM for [a| < k — 2 in boundary normal coordinates. 

4. Local interior rigidity; Proof of Theorem [1] 
Given a symmetric 2-tensor / = , the divergence of / is an 1-tensor 5f defined by 

[Sf]i = !l' h V,f„ 

in any local coordinates, where are the covariant derivatives of the tensor /. Given an 1-tensor 
(a vector field or an 1-form that we identify through the metric) v, we denote by dv the 2-tensor 
called symmetric differential of v. 

[dv]ij = - (ViVj + VjVi) . 
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Operators d and —5 are formally adjoint to each other in L 2 (M). It is easy to see that for each 
smooth v with v = on dM, we have I(dv)('y) = for any geodesic 7 with endpoints on dM. This 
follows from the identity 

(13) ^M7(t)),7(*)> = <*(7(t)),7 a (t)>- 

It is known (see [Shl| and (|14p below) that for g smooth enough, each symmetric tensor / € 
L 2 (M) admits unique orthogonal decomposition / = f s + dv into a solenoidal tensor Sf := f s 
and a potential tensor T 3 / := dv, such that both terms are in L 2 (M), / s is solenoidal, i.e., 5f s = 
in M, and u € Hq(M) (i.e., t> = on dM). In order to construct this decomposition, introduce 
the operator A s = 5d acting on vector fields. This operator is elliptic in M, the Dirichlet problem 
satisfies the Lopatinskii condition, and has a trivial kernel and cokernel. Denote by A S D the Dirichlet 
realization of A s in M. Then 

(14) v = (A s D )- 1 5f, f = f-d(A s D y l 5f. 
Therefore, we have 

V = d(A s D y 1 d, S = ld-V, 

and for any g E C (M), the maps 

(15) (Afj) -1 : H~\M) -»■ H^(M), T,S : L 2 {M) — > L 2 (M) 

are bounded and depend continuously on g, see |SU31 Lemma 1] that easily generalizes for manifolds. 
This admits the following easy generalization: for s = 0, 1, . . . , the resolvent above also continuously 
maps H^ 1 into H s+l n Hq, similarly, V and S are bounded in H s , if g G C k , k S> 2 (depending 
on s). Moreover those operators depend continuously on g. 

Notice that even when / is smooth and / = on dM, then f s does not need to vanish on dM. 
In particular, f s , extended as to M, may not be solenoidal anymore. To stress on the dependence 
on the manifold, when needed, we will use the notation vm and ff^ as well. 

Operators S and V are orthogonal projectors. The problem about the s-injectivity of I, restricted 
to a subset of geodesies, then can be posed as follows: if If = on those geodesies, show that 
f s = 0, in other words, show injectivity on the subspace SI? of solenoidal tensors. 

We start with the proof of Theorem [H We will split the proof into several steps. 

4.1. Choosing a suitable metric isometric to g. Any two metrics such that one of them is a 
pull-back of the other under a diffeomorphism fixing the boundary pointwise, will be called below 
isometric. Such a diffeomorphism is necessarily C k+1 if the metrics are C k , see e.g., |SU3j . and the 
norm of its derivatives are controlled by those of the two metrics, see [STJ3], Lemma 6]. 

We first find a metric isometric to g, that we denote by g again, so that the boundary normal 
coordinates related to g and g coincide in some neighborhood of the boundary. 

By |CDSl Theorem 2.1], if e <C 1, there is h € C k ~ l isometric to g so that h is solenoidal w.r.t. 
g. Moreover, 

(16) \\h-g\\ ck -t<Ce, 

and there is no need to replace k by k — 1 if we work in the C k ' a spaces. By a standard argument, 
by a diffeomorphism that identifies normal coordinates near dM for h and g, and is identity away 
from some neighborhood of the boundary, we find a third g\ isometric to h (and therefore to g), so 
that gi = g near dM, and g\ = h away from some neighborhood of dM (and there is a region that 
gi is neither). Then g\ — h is as small as g — h, more precisely, 

(17) ||<?i - h\\ c i-3 < C\\g — h\\ c i-i, l<k. 
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This follows from the fact that gi 
and 



*h, with a diffeomorphism <f> that is identity on the boundary, 
- Id|| c ;-2 < C\\g 



h\ 



f 



f = h-g, f = gi 
0. Estimate (|17p implies 

llZ-ZllcJ-s^cil/lb-i, 



C' k 



C i-x. 



(J- 



VZ < k. 



< Ce. 



(18) 
Set 
(19) 

We aim to show that / 
(20) 

By ©, CED and {2D]), 

(21) ||/||c*-i < Ce, 
By Theorem [3l and the remark after it, 

(22) d a f = on dM for \a\ < k — 5. 

We have now two isometric copies of g: the first one is h that has the advantage of being solenoidal 
w.r.t. g\ and the second one g\ that has the same jet as g on dM. We need both properties below 
to show that g = h, i.e., / = (or g = g%, i.e., / = 0) but so far we cannot prove that h = g\. The 
next proposition shows that h and g\ are equal up to 0(||/|| 2 ) = 0(||/|| 2 )- 

Proposition 1. Let g and g be in C k , k > 2 and isometric, i.e., 

g = i>*g 

for some diffeomorphism ip fixing dM pointwise. Set f = g — g. Then there exists v vanishing on 
dM, so that 

f = 2dv + f 2 , 

and for g belonging to any bounded set U in C k , there exists C(U) > 0, such that 

IIMIc*-' < c(u)\\^ - id||^ fc _ x , |M| c *- a < c{u)\\^ - id|| c *_i. 

Proof. Extend g to M in such a way that the C k norm of the extension is bounded by CH^Hc^M)- 
Set v(x) = exp~ 1 (tp(x)) that is a well defined vector field in C 1 (M) if ip is close enough to identity in 
C 1 (it is enough to prove the theorem in this case only), and u = 0on dM. Set ip r (x) = ex.p x (rv(x)), 
< t < 1. Let g T = ip*g. Under the smallness condition above, v is small enough in C , and 
therefore ip T is close enough to identity in the C l {M) norm. Therefore, tp T : M — > ip T (M) C M is 
a diffeomorphism. Next, tp T fixes dM pointwise, therefore, ip T {M) = M. 
The Taylor formula implies 



d 

9 + Tr 



where 



T=0 



1 



g T + h = g + 2dv + h, 



\h\ < — max 

2 re[o,i] 



d 2 g T 



dr 2 



and 2dv is the linearization of g T at r = 0, see [Shi] . To estimate h, write 

dij) k difj l T 



9ij = 9ki ° ipr 

and differentiate twice w.r.t. r. Notice that 

d 2 ip. 



dr 2 



<C\\v\\ 2 L c 



dx l dxi ' 
d 2 Vip T 



dr 2 



< C\\v\\ 2 cl 
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This yields the stated estimate for fi for k = 2. The estimates for k > 2 go along similar lines by 
expressing the remainder h in its Lagrange form, and estimating the derivatives of h. □ 

We apply Proposition [1] to h and g\ to get by (fT8|) . 

(23) f = f + 2dv + f 2 , ||/ 2 || C «-»<C||/||2 7l _ 1> VI <k. 

In other words, f s = f up to 0(||/|[ 2 ). 

Next, with g extended as above, we extend g\ so that gi = g outside M. Then g G C fc and 
5i G C fe " 5 by (E2D. 

4.2. Reparametrizing the scattering relation. We proceed with some preliminary work that 
would allow us to apply |SU4[ Theorem 2]. Assume first that the underlying metric is fixed to 
be go. Let (M, go) be a C k extension as above. In [SU4] . the geodesies are extended to M \ M, 
parametrized by initial points, and corresponding directions on a finite collection {H m } of smooth 
connected hypersurfaces in M, having additional properties as explained below. Given two complete 
T>' <s T>, we will construct such a family issued from a set T>" with V <e T>" <e T>, that is also 
complete. 

For any (zq,(o) G T*M, including the case where zq G dM, there is a maximal geodesic 70 
through zq normal to Co> satisfying the conditions of Definition [2j Let us assume that 70 is 
parametrized by t G ±P > 0, and 7(0) = z, 70(0) = £o> with conormal to £o- Since 70 

is maximal in M, by Lemma [TJ for any 5± > there exists 5 G (0, <5i) so that the extension 70 of 7 
to M corresponding to t G [l~ — 5, l + + 5] is well defined and has endpoints in M mt \ M. By ([6]), 

{( 7o (t),7o(0); r <t < Z + } n5(9M) c v. 

The extension (70,70) may have additional point on dM corresponding to l~ — 5 < t < l~, and 
l + < t < l + +5. However, we choose 5i <C 1 so that they still belong to T>, and this is possible to do 
because V is open. Now, by LemmaEl any geodesic that is obtained by a small enough perturbation 
(z, £) of the initial conditions (zq, Co) of 70 at t = 0, with the same interval t G [l~ — 5, l + + 5], will 
satisfy condition ([9]). Condition ([7]) will also be satisfied by a perturbation argument, if 5 is small 
enough as well. Now we can perturb g$ in the C 2 (M) topology to ensure the same property. Note 
that if \\g — go\\c 2 (M) — e > one can choose an extension g of g to M so that \\g — go\\c 2 (M) — C £ 
with 50 a fixed extension of 50 as above. 
To summarize, we proved the following. 

Lemma 4. Under the conditions of Theorem^ for any (zo^Co) £ S*M, there exists a geodesic in 
the metric go through zq normal to C,q with endpoints in M mt \ M so that conditions (OJ) are 
satisfied. 

Moreover, if that geodesic has initial conditions (zo>£o) at t = to for some to G [0,1], and an 
interval of definition < t < I, properties (?]) and (0) remain true under small enough perturbations 
of (z ,Co), and g in C 2 {M). 

Let us assume now that the underlying metric is g as in the theorem, with e<1. Since T>' is 
compact, there are finitely many geodesies 

with the following properties. If 7 m (0) = £ m G S Zm M, then for any m there exists neighborhoods 
U' m (s U m of (z m ,£ m ) in SM, such that if T m , T' m is the set of geodesic with initial conditions in 
U m , respectively U' m , and the same interval of definition as r y m , then for any (x, Q G T*M there is 
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a geodesic 7 G UT^ so that ((7J), © are satisfied with V replaced by V; and all geodesies in ur m 
satisfy ((JJ), ([9|) as well. Moreover, 

(24) uCD^fF), 

where is regarded as a point set in the phase space SM consisting of the points on all integral 
curves. 

We will parametrize T m , T' m , with initial points outside M. Choose a family of finitely many 
small enough smooth hypersurfaces {H m } in M mt \M, each one transversal to 7 m . Without loss of 
generality, we can assume that all geodesies in T m can be extended in M mt \M so that they intersect 
H m once and are transversal to H m as well. We still denote the set of the extended geodesies by T m 
and T' m , respectively. Let Ti m be the open subset of € SM; x G H m , £ g" T x H m } formed by 

those (x, £) that coincide with the left endpoint and the corresponding direction of some geodesic 
in T m . Then their endpoints belong to M mt \ M again, and their length is a smooth function 
lm(%,Q > (actually, we may even assume that l m is constant, if U m is small enough). 

We have 

(25) T m = T(H m ) = {7«,€(t); < t < l m (x,£), (x,0 G H m } , 

where 7* £, as usual, is the geodesic with initial conditions (x,£). We define Ti' m in a similar way, 
related to U' m . We consider also the geodesies in the metric g defined as in (f25l) . Then ((2J), ([9]) hold 
for those geodesies, too, provided that e <C 1. 

Combining the arguments above with Lemma EJ we have the following. 

Proposition 2. Let g and g\ be as in the end of Section \4-l\ Then, if e 1, 

<S> l ™( x >t\x,0 = ¥™ ix 't\x,0, V(x,0 G Vm, 

where is related to g\ . Moreover, ur m satisfies ^j, (GJ), a^rf is complete in the sense that 
N*(\JT' m ) D T*M, and satisfies {2J$ ; similarly L)T m and UT' m have the same properties. 

In other words, informally speaking, we pushed the boundary, where the scattering relation is 
defined, to a collection of hypersurfaces outside M, so that the corresponding geodesies are always 
transversal to them, and the endpoints are away from dM. 

Proof of Proposition^ The proof is straightforward, if the geodesic issued from (x,£), for < I < 
l m (x,£), always intersects dM transversally. Observe that k± = fc± because g and g have the same 
normals on dM, therefore = implies S = E. The points (x,£) G Hm where this transversality 
does not hold is a closed set of measure zero by Lemma 13 and for such points, one can approximate 
with points outside this set, and to use the continuity of <5*. □ 

4.3. Linearization near g. Now we are in the situation of [SU4J, sec Theorem 2 there. Choose 
smooth functions a m supported in TC m and equal to 1 on H' m . Set a = {a m } and I am = a m I, more 
precisely, 

la m f(z,0= a m(x,0 {f(j z £),J%t)dt, (z,^)£ Hm- 

JO 

Also set 

la = {Ia m } , N am = I* m I am , N a = N am , 

where the adjoint is taken w.r.t. the measure \{v, £)|dE2n_2, where dE2 n -2 is the induced measure 
on dSM by the volume form, and v is a unit normal to H m . 
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In |SU4|, Theorem 2] we showed that if for a fixed g, I a is s-injective, then we have the following 
a priori estimate 

(26) \\f s \\L^M)<c\\N a f\\ mm , 

with a suitable norm || • ||#2(m) so that H 2 (M) C H 2 (M) C il^M). Moreover, (f26|) remains true 

under small C k perturbations of g with a constant C that can be choose uniformly. Note that (T) 
is not needed in [SU41 Theorem 2]. 

Fix m and £ ft m . Let {7a,,e(*), < i < Z m (a:,£)}, and {^(t), < t < / m (x,£)} be the 

geodesic issued from (x,£) related to g and gi, respectively. By Proposition [21 their endpoints and 
directions coincide. We reparametrize 7^, so that t £ [0, 1]; then they have the same speeds 

l7«,f(t)| = = l m(x,0 ■ 

Define the following variation of 7x,£> where exp is related to g: 

(27) c T {t) = exp 7 ^ (t) (Tv(t)) , v(t) = exp-^ {f) (%^(t)) , 
where < t < 1, < r < 1. Then c$ = 7^, c\ = %£■ Set g T = g + r(g\ — g). Let 

(28) E(t) = [ {c T (t),c T (t)) 9T dt, 

Jo 

where, in local coordinates, {c T (t),c T (t)) g T = gJ^{c T )ci r 8 r . Apply Taylor's formula 

E(l) = E(0) + E'(0) + I (1 - t)E"(t) dr 

to get 

(29) E'(0) = - f (1 - t)£"(t) dr 
because £(0) = = l 2 m {x,^). Write 



*1>(t,s)= {c T (t),c T (t)) g sdt= 4(c r )44dt, 

JO JO 

where the second integrand is written in local coordinates. Then E(t) = iP(t,t). For E' we get 
(30) E'(t) = Mt,t) + Mt,t). 

Since Co = 7^ is a critical curve for the energy functional, we get V>s(0, 0) = 0, therefore, 

E'(0)= / 1 (/,7^)di, 



(l-T)E"(r)dT. 



recall flU]). Together with ([29]) this yields 

(31) //(7a,.*) = - 

To estimate the r.h.s. above, note that 

(32) E"(t) = ipfT (r, t) + 2^ rs (r, r) 
because ip ss = 0. Note that 



(33) 



dc T (t) 



+ 



dCr(t) 



<C(\v(t)\ + \v(t)\), 
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(34) 



d 2 c T (t) 



3r 2 



+ 



d 2 t T {t) 



dr 2 



<C(\v(t)\ + \v(t)\) 2 . 



We have that | exp" 1 y\ < C\x — y\ for \x — y\ <C 1, where the norm is in any fixed coordinate chart, 
and C > depends on an upper bound of g in C k , > 2. All constants below will have the same 
property. This and (|27p imply 

(35) \v(t)\<C\%^t)- lx ^t)\<C'\\f\\ cl . 

Since in fixed coordinates, D x exp" 1 y + D y exp" 1 y = when x = y, we have 

\D x exp x 1 y + D y ex.p~ 1 y\ < C\x-y\. 
This allows us to estimate v(t), see (|27|) . as follows 

(36) \v(t)\ <c(\%^t)- lx4 (t)\ + \%^t)-% A (t)\) <C'\\f\\cx. 

By iU), (EH), (J33D, dMD, §5}, dSHD, 

|//(7^)I<C'||/||^. 

This is the same estimate that was used in the linearization argument in [StJ2~[ ISU3| and goes back 
to jE], but now proven in this more general situation. Therefore, 

(37) <C||/||^. 
That implies the same for HiVa-ZH^oo, see e.g., |SU3| . therefore, 

(38) \\N a .f\\ L o <C\\f\\ 2 cl , Vj. 

4.4. End of the proof of Theorem [H We will use interpolation to estimate \\^ajf\\H 2 (M) 
through some power of \\N a . /Hi 00 - Since N a . is a ^DO of order —1 and /, extended as outside 
M, is smooth enough if k 3> 2 by ([22]) . we get 

(39) W^M^m) < ll^/f ff2( M) ^ C\\f\\ mM) \\N a J\\ L2{Ay 
Combine d3HJ) and ((39]) to get 

(40) W N JWh(M) ^ c|l/lb3||/ll 2 cl < c'\\ff c5 , 

by m- 

Since I go ,v is s-injective, so is N a , related to go, by the support properties of a. Now, since N a 
by (JH]) is s-injective for g = go, we get from [SU41 Theorem 2] that N a (the one related to g) is 
s-injective as well provided that 

( 41 ) \\9 -9o\\c{M) ^ £ 

with some k 3> 1 and £o > 0. Moreover, (f26|) is true. Assume now that both (|4ip and (jHJ) 
are satisfied. The geodesies issued from suppa form a complete set by the second statement of 
Proposition El therefore, by (jUJ) and ([26]), 



(42) Wf'Wvm < c\\nJ\\ &2 < c'\\ffjl 

By mJ s = f + fl, therefore, 

||/ s ||l2 (M ) > ||/|| L 2( M) - ||/|||l2( M) > ||/||l2 (M ) - C\\f\\ 2 C 2. 
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Together with (|4"2"j) . this yields 

II/IMm)<c(||/|| 2 C2 + 11/11^) <C'\\f\\% 

because the C 5 norm of / is uniformly bounded when e < 1. We can use now interpolation estimates 
in C k , see |Tri| . and Sobolev embedding estimates to get ||/||e?s < C||/||^ < C"||/||^. n/2+1 with any 
H € (0, 1) as long as k = S> 1 in ([2T|) , Next, interpolation estimates in Sobolev spaces imply 
H/ll^-n/2+i < C||/lli2) so in the end, we get 

3ti 2 /2 



L2(M) < C \\f\\ L 2 {M y 

Choose 2/3 < /i 2 < 1 with a corresponding k S> 2, to get ||/||l 2 (m) > 1/C if / ^ 0. Choose e< 1 
to get a contradiction with (|2ip . This proves Theorem [2] for e replaced by min(eo,e)- 

Now, / = implies h = g, therefore, g and g are isometric. 

This concludes the proof of Theorem [TJ 



5. Proof of Theorem [2] and Corollary Q] 

Proof of Theorem^ As we mentioned in the Introduction, Theorem [2] is a reformulation of [SU4, 
Theorem 3], as we show below. Notice first, that the s-injectivity of I g .x>" f° r some T>" C V implies 
s-injectivity of I g T>'- One can always assume that M is equipped with a finite analytic atlas. Note 
that the assumption that V is open implies that the corresponding set of geodesic is open in the 
sense of [SU4] by Lemma (2) Let T>" (s T>' be such that T>" is still complete and open. It can be 
constructed as in Section 14.21 As in Section 14.21 again, with T>" and T>' playing the roles of V and 
T> there, respectively, we construct N a such that all geodesies through suppa cover kZ 1 {T>") and 
are contained in the interior of kC 1 (£>'). Then N a is s-injective for each analytic go € Q by [SU4, 
Theorem 1]. It is still s-injective under a small enough C k , k » 2, perturbation of g G C k {M) by 
[STJl Theorem 2]. Note that [SUil Theorem 2] requires that the perturbation must be considered 
in C k (M) but one can use extensions of g near a fixed go so that their norms in C k (M) are bounded 
by llsilofc(M) with a fixed C. Using the fact that V" <g V , LemmaO and the support properties of 
a, we deduce that / 5i x>' is s-injective for g close enough in C k {M) to a fixed analytic go S Q- We 
can therefore build Q s as a small enough neighborhood of the analytic metrics in Q, that form a 
dense set. 

This completes the proof of Theorem [21 □ 

Corollary [1] is an immediate consequence of Theorem [1] and Theorem [2l 

Proof of Remarkm Condition (T) is not needed for [SU4, Theorem 2], see also (|26p but it is used 
in the proof of Theorems 2 (s-injectivity for analytic metrics) and Theorem 3 (generic s-injectivity) 
there. Assume now that (T) in Theorem [2] of this paper is replaced by the assumption that (M,g) 
is simple as in the remark. In the proof above, in this situation, we need to show that N a is 
s-injective for a dense set of metrics in Q. Fix go € Q. It can be extended as <7o to M so that (M,g) 
is simple. Given e > we can find a real analytic g so that \\g — 9o\\ck(M) — £ f° r ^ fixed. Then 

the ray transform / related to g over all maximal geodesies in M is s-injective, see |SU3|. [SU4] , Let 
now Ng ta f = in M, where g = §\m, and the subscript g in N g ^ a indicates that this is the normal 
operator related to g. Then we get after integration by parts that I g f{^) = for all maximal 
geodesies in M. Let / be the extension of / as zero outside M. Then I g f(j) = for all maximal 
geodesies in M. Therefore, / = dv in M, v E Hq(M). Since / = in M \ M, on can see that the 



LENS RIGIDITY 



17 



same holds for v as well, see also [Sh2| . Therefore, / is potential, thus I g is s-injective, and so is 
Ng, a , see jSU4[ Lemma 2]. □ 
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